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Heat transport in low-dimensional systems has attracted enormous attention from both theo- 
retical and experimental aspects due to its significance to the perception of fundamental energy 
transport theory and its potential applications in the emerging filed of phononics: manipulating 
heat flow with electronic anologs. We consider the heat conduction of one-dimensional nonlinear 
lattice models. The energy carriers responsible for the heat transport have been identified as the 
renormalized phonons. Within the framework of renormalized phonons, a phenomenological theory, 
effective phonon theory, has been developed to explain the heat transport in general one-dimensional 
nonlinear lattices. With the help of numerical simulations, it has been verified that this effective 
phonon theory is able to predict the scaling exponents of temperature-dependent thermal conduc- 
tivities quantitatively and consistently. 

PACS numbers: 05.45.-a,05.60.-k 



I. INTRODUCTION 



The Fourier's heat conduction law states that the heat current flowing through a system is proportional to the 
temperature gradient imposed on the two ends of the system: j ~ — kVT, where j is the heat current, VT is the 
temperature gradient and the proportionality k is defined as the thermal conductivity. A diffusive or normal heat 
conduction requires that k is independent to the system length L in the thermodynamical limit L — > oo. However, 
the discovery of anomalous heat conduction [l| that k cx L" with < a < 1 for one-dimensional Fcrmi-Pasta-Ulam 
/3 (FPU-/3) lattice has casted doubt on the validity of Fourier's heat conduction law on low-dimensional systems and 
stimulated intensive studies on this issue from both theoretical [H-Q and experimental 0] attempts. On the other 
hand, it has been demonstrated that the nonlinearity can be utilized to design novel nanoscale solid-state thermal 
devices such as thermal diodes [8l4ll| . thermal transistors [l^, thermal logic gates [l^ and thermal memories [13, [iBl 
which gives birth to the innovating field of phononics: manipulating/controlling heat fiow and processing information 
with phonons [l^. The distinctive and unique transport property of low-dimensional system has posted great challenge 
to the complete microscopic transport theory. Therefore, any theoretical attempt towards a thorough understanding 
of the heat transport in general one-dimensional nonlinear lattice systems is timely and highly desirable. 

In order to reveal the physical mechanism underlying the heat transport in one-dimcnsional nonlinear lattices, the 
energy carriers responsible for heat transport must be identified on the first place. Although the linear Harmonic 
lattice can only sustain the vibrations of phonons, it has been known that there are more than one typ e of excitation 
modes in nonlinear lattices, i.e. the renormalized phonons [T7l - |2^ . solitons [25l - [28j and breathers [29l l30l| . In particular, 
the stable solitons which almost do not interact with each other has been argued to be the origin of the anomalous 
heat conduction found in FPU-/3 lattice [sH - fssj . The numerical calculations for the sound velocities of energy carriers 
have also been found to follow the prediction of soliton velocities [12] ■ However, these calculations cannot exclude the 
possibility of renormalized phonons as en ergy carriers since the prediction of sound velocities of renormalized phonons 
is not far from the prediction of solitons [23|. It will be necessary to determine which excitation mode is the energy 
carrier for nonlinear lattices and more accurate calculations need to be performed to uncover this confusion [sj- 

For the past decades, enormous efforts have been focused on the study of normal or anomalous heat conduction, 
i.e. size-dependent thermal conductivities of one-dimensional nonlinear lattices [l|-[^ . On the contrary, there are only 
a few works dealing with the temperature-dependent behaviors of heat conduction [H, [l^, [H, [35l - [4]| which should 
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be more relevant to experimental investigations. The temperature dependence of heat conduction is caused by the 
nonlinearity as the linear Harmonic lattice cannot display temperature modulated response. Furthermore, the novel 
thermal devices such as the thermal diodes are rooted on the fact of temperature or nonlinearity modulated phonon 
spectrum. The understanding of the temperature-dependent behavior of heat conduction for nonlinear lattices is 
thus of both theoretical and experimental importance. A transport theory based on the correct energy carriers will 
definitely help us to unravel the underlying physical mechanism for the heat transport in low-dimensional nonlinear 
lattices. 

The paper is organized as the follows: in Sec. II the general one-dimensional nonlinear lattice models will be 
introduced and the the concept of renormalized phonons will be discussed. The sound velocities of energy carriers will 
be analyzed and determined through detailed numerical calculations. Sec. Ill will then present the phenomenological 
theory of effective phonon theory in the framework of renormalized phonons. The predictions for temperature- 
dependent thermal conductivities from effective phonon theory will be compared with Non-Equilibrium Molecular 
Dynamics (NEMD) simulation results. We will give conclusions and summaries in Sec. IV. 



II. NONLINEAR LATTICES AND RENORMALIZED PHONONS 



The one-dimensional nonlinear lattice Hamiltonian has a general form: 
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where pi and Xi denote the momentum and deviation from equilibrium position for the i-th atom, respectively. For 
simplicity, the periodic boundary condition xi = x^+i is often applied where is the number of atoms or the size of 
lattice in dimensionless unit since L = Na where a is the lattice constant and can be set as unity [l^. The inter-atom 
potential V{xi^i^Xi) and the on-site potential U{xi) also assume the following general expressions: 
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For simplicity, the coefficients Vg and Us only take two values and 1. The celebrated FPU-/3 lattice 
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has two nonvanishing coefficients V2 = i'4 = 1 while the most studied on-site 0** lattice 
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(4) 



has two nonvanishing coefficients with V2 = un ~ \. 

For the Harmonic lattice, it is well-known that the Hamiltonian can be decomposed into the energy sum of normal 
modes (phonons) under the transformation qk — X^i^i ^kiXi and = ^kiPi where the matrix clement Ski — 

-^{sm2Trki/N + cos2TTki/N). The Hamiltonian in position space can thus be transformed into normal mode space 



as 
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where is the phonon frequency. The resulted total Hamiltonian in phonon space is the energy summation of every 
independent phonon modes. Therefore, the Harmonic lattice has no phonon-phonon interaction and is called linear 
lattice in this sense. 

As for the nonlinear lattice, there does not exist any transformation which can transfer the lattice Hamiltonian into 
the sum of independent phonon modes. However, the nonlinear lattice Hamiltonian can be expressed approximately 
as the energy summation of phonon modes with renormalized frequencies due to the nonlinearity. 
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The connection between the renormaUzed phonon frequency and the nonhnearity originates from the generahzed 
equipartition theorem: 

dH \ 



where ks is the Boltzmann constant and (•) denotes ensemble average. For Harmonic lattice, it is easy to see from 
Eq. (O that 

kBT - ul (ql) (7) 

which is the familiar result that every degree of freedom (phonon mode) shares the same amount of energy. For 
nonlinear lattice, a similar relation holds for the renormalized phonon modes under the mean-field approximation [2^ : 

kBT « UjI (ql) (8) 

where the renormalized phonon frequencies can be defined as 
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and the renormalization coefficients a and 7 encode the information of inter-atom potential and on-site potential, 
respectively. For the simple Harmonic lattice with i>2 = 1, the renormalization coefficient a = \ and 7 = 0, the 
renormalized phonon frequencies Wfc = oj^. thus recovers the original phonon dispersion relation. 

The nonlinear lattices can be classified into two categories: with or without on-site potential. For lattices without 
on-site potential as U{xi) = 0, we have 7 = from Eq. ([9]). Therefore, the branch of renormalized phonon modes is 
acoustic-like with 

Wfc = y/a ■ Wfe (10) 

For lattices with on-site potential as U{xi) 7^ 0, the renormalization coefficient 7 7^ and the resulted phonon branch 
is optic-like. 

In order to understand the physical mechanism of the heat transport, wc need first to determine the energy carriers 
in nonlinear lattices. Besides the renormalized phonons, other excitation modes such as solitons and breathers existing 
in nonlinear lattices are also possible candidates for energy carriers. One of the best way to identify the energy carriers 
is to calculate the temperature-dependent energy transport speeds and compare them with the theoretical predictions 
for the sound velocities of different excitation modes. 

The sound velocity of renormalized phonons can be derived from Cs = ^^\q=o where we have applied the substitution 
q = 2Trk/N and = 2 smq/2 in the continuous limit. For the well-known FPU-/? lattice with V2 = I'i ~ 1, the sound 
velocity of renormalized phonons can be obtained from Eq. ([5]) , ([9]) and ([T0|) : 




'da;-a;4e-(^'/2+^V4)/T 

(11) 



where dimensionless unit has been used and the Boltzmann constant ks has been set as unity. The sound velocity 
depends on the temperature, or equivalently the strength of nonlinearity. In the low temperature limit T — > 0, — > 1 
as expected for Harmonic lattice. In the high temperature region T » 1, the sound velocity Cs w 1.22T-'^/* exhibits 
strong nonlinear dependence. On the other hand, the soliton theory predicts a temperature dependent sound velocity 



clyM^ = 7^T (12) 

where 77 is a fitting parameter. By tuning 77 = 2.215, the sound velocity from soliton theory of Eq. (|12p coincides 
with the prediction from renormalized phonons of Eq. (|lip in low and high temperature region. However, these 
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FIG. 1: (color online). Sound velocity Cs as the function of temperature T for FPU-/3 lattice. The circles are the numerical 
data. The dotted line is the theoretical prediction from renormalized phonons from Eq. The dashed line is the theoretical 

prediction for solitons derived from the formula of Eq. (|12|l with rj — 2.215. The cross symbol is the numerical data for T = 0.5 
obtained in Ref. Adapted from Ref. 



two predictions deviate from each other in the intermediate temperature region and can only be verified by accurate 
numerical calculations [42j. 

In Fig. [TJ the numerically calculated sound velocity Cs for FPU-/? lattice has been compared virith the prediction from 
the theory of renormalized phonons and solitons [3^, respectively. It can be clearly seen that the soliton prediction 
deviates from numerical data in the intermediate temperature region while the prediction of renormalized phonons 
matches perfectly with the numerical data in all temperature regions. This is a clear evidence that the energy carriers 
in the FPU-/? lattice are the renormalized phonons rather than the previously believed solitons. 

To demonstrate the validity and consistency of the renormalized phonon formulation, the other three iJ„ models 
with n = 3, 4, 5 have been considered with the following Hamiltonian: 



N 



(13) 



The Hi model is just the high temperature limit of FPU-/? lattice. In the framework of renormalized phonon formu- 
lation, the sound velocities of Hn model can be expressed in a compact form: 



x"-e 



(14) 



In Fig. [21 the numerically calculated sound velocities are compared with the predictions for renormalized phonons 
for Hn lattices with n = 3,4, and 5. Again, the quantitative agreements have been found for all three models which 
we considered. 

The above examples are for the lattices without on-site potential. For the lattices with on-site potential such as 0'* 
lattice, there are numerical results suggesting that the energy carriers are also renormalized phonons [43j. 
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III. EFFECTIVE PHONON THEORY 



So far as wc know, there is no microscopic transport theory which is able to explain the heat transport in nonlinear 
lattices with strong nonlinearity. The perturbation theory will become invalid as the nonlinear term cannot be 
ignored anymore in this case. For example, no existing theory can even predict the temperature dependence of 
thermal conductivities for general nonlinear lattices. A phenomenological theory possessing such kind of prediction 
power is thus highly desirable. 

The thermal conductivity for solids is usually described by the Deybe formula: 



c 



- dqv'^r, (15) 



where c denotes the specific heat, Vq denotes the group velocity and Tq denotes the relaxation time for phonon mode 
q. The key issue for the Deybe formula is then to determine the relaxation time Tq. 

For nonlinear lattices, we have shown the evidences that the energy carriers should be the renormalized phonons. 
Therefore, the transport theory should be constructed with the basis of renormalized phonons. In consideration of 
the nonlinearity and anomalous size-dependence of thermal conductivity, we propose the modified Deybe formula in 
the framework of renormalized phonons: 



'^=7r dqv^TqPiq) 

Jo 



(16) 



where P{q) is a weight factor with normalization condition J^^ dqP{q) = 1, is the group velocity of renormal- 
ized phonons as Vq — du)q/dq, and Tq is the relaxation time for renormalized phonons which assumes the following 
proportionality: 

r.ocl.?^ (17) 



where coq is the frequency for renormalized phonons and the dimensionless parameter e represents the strength of 
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nonlinearity defined as: 

, \{En)\ 



{El + E„) 



0<e<l (18) 



where Ei denotes the hnear potential energy and En denotes the nonlinear potential energy. The parameter e describes 
the ratio of nonlinear potential energy to the total potential energy and is independent of mode q. 

The introduction of the weight factor P{q) is to be consistent with the anomalous size-dependent thermal con- 
ductivities for lattices without on-site potential such as the FPU-/? lattice. If we assume P{q) oc 1/q^ in the long 
wave-length limit g — )- 0, it is straight forward to obtain the size-dependent k as 

K (X dqVgTqP{q) cx / dq — • — cx q~^ cx (19) 

Jo+ Jo+ Q 

the lower integration limit O'^ describes the fact that the lattice size is always finite with length N and oc g cx 1/A^ 
in the long wave-length limit g — due to the finite size effect. The exact value of exponent 6 is still unresolved and is 
argued to be 2/5 from mode-coupling theory or 1/3 from hydrodynamic theory Q. However, there are one exception 
that the momentum conserved one-dimensional coupled rotator model without on-site potential does exhibit normal 
heat conduction behavior [3 |45|. Most recently, some numerical works even suggest that the momentum conserved 
one-dimensional lattice models with asymmetric interactions might also follow the Fourier's normal heat conduction 

lawjii,!!!. 

It should be noticed that the introduction of P{q) simultaneously explains the normal heat conduction for lattices 
with on-site potential. The non-zero renormalization coefficient 7 shifts the renormalized phonon frequency from zero 
to ujq = y/j in the long wave-length limit g — >■ 0. The integral in last equation becomes 

2ir 1 1 1'^'" 



K (X / dq- ^ < — / dq^ (20) 

Jo+ T Vl Jo+ T 

One should notice that S < 1 due to the constriction of normalization condition for P{q) as ^ J^^ dqP{q) = 1. The 
integral is non-divergent and the thermal conductivities for lattices with on-site potential is finite and obeys Fourier's 
heat conduction law which is consistent with previous results jsil [35ll48j. 

In this paper, we will only focus on the temperature dependent thermal conductivities for nonlinear lattices without 
on-site potential. The renormalization coefficient 7 = and the renormalization phonon frequency Coq = y/aujq where 
temperature dependent part a and mode dependent part ujq are separate. From Eq. (|16l) and (|17p . the temperature 
dependence of thermal conductivity can be fully described by the simple formulation 

^{T) oc ^ cx ^ (21) 

where the renormalization coefficient a (or the sound velocity Cs) and nonlinearity strength e are fully determined by 
the lattice Hamiltonian. One should also notice that these two parameters only depend on temperature and are mode 
independent. In another word, the size- and temperature-dependences of the thermal conductivities are separate, 
which enables us to discuss the temperature-dependence of the thermal conductivities alone without worrying about 
the annoying size-dependence. This effect has been verified by our numerical simulations (not shown here). 

For the FPU-/? lattice with Hamiltonian of Eq. the temperature dependence of Cg and e exhibits scaling 

relations at both low and high temperature limits. As we have discussed above, the sound velocity Cs shows the 
following temperature dependence: 

Cs CX const, T « 1; Cs cx T^/'', T » 1 (22) 
From definition of Eq. (fT8|). the nonlinearity strength e for FPU-/? lattice can be expressed as 

(E,(x,+i-xOV4) 



(23) 



At low temperature limit, the quartic potential term is a small term comparing to the quadratic potential term. 
The whole system can be described by a harmonic lattice plus a small perturbation. From equipartition theorem for 
harmonic lattice, we thus have {J^ii^i+i ~ ^iY) ~ NksT for the FPU-/? lattice at low temperature limit. This also 



7 




FIG. 3: (color online). Thermal conductivity k as the function of temperature T for FPU-/3 lattice. The lattice size is set as 
N = 1000. The circles are the numerical data from NEMD simulations. 



gives rise to the scaling relation for the quartic potential term, i.e. {^^ixi+i — x;)^) oc T^. Thus the nonlinearity 
strength e follows 

On the other hand, the quadratic potential term is a small term at high temperature limit and the nonlinearity 
strength e approaches to the upper limit 1. The entire temperature dependence of e can be summarized as: 

e cx T, T « 1; e oc const, T » 1 (25) 

From Eq. (f2T|) . ([22|) and ([25|) . we can derive the temperature dependence of k for FPU-/3 lattice at both low and 
high temperature: 

K(T)oci T<<1; K(^)ocTl/^^»l (26) 

In Fig. in the numerical results of thermal conductivities k calculated by Non-Equilibrium Molecular Dynamics 
(NEMD) simulations are plotted for FPU-/3 lattices. It can be seen that the scaling of the thermal conductivities 
follows the relations k oc l/T at low temperature limit and n oc T^/^ at high temperature limit, respectively. The 
temperature dependence of k for FPU-^ lattice can thus be consistently explained by our effective phonon theory at 
both low and high temperature limits. 

To further verify the predictions from effective phonon theory, we also consider the i7„ lattices with Hamiltonian 
of Eq. (fT3)) . The fact that the potential energy of these iJ„ lattices is fully nonlinear gives rise to the special property 
for the nonlinearity strength as e = 1. From Eq. (|14p . the temperature dependence of the sound velocity Cg follows 
Cs oc yi/2~i/'i_ The thermal conductivities n for lattices can be derived as: 

KocCsOcTs-^ (27) 

Therefore, from effective phonon theory, k oc 7^1/6^7^1/4 g^^j jnS/io j^j. jj^^ ^^^^ jj^ lattices, respectively. From 
Fig. m the numerical results of thermal conductivities for all these iJ„ lattices are plotted and compared with the 
theoretical predictions from effective phonon theory. Perfect agreements between theory and numerical results have 
been found for all the _ff„ lattices at temperature regions over three orders of magnitudes. 
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FIG. 4: (color online). Thermal conductivity n as the function of temperature T for Hn lattices with n = 3,4 and 5. The 
lattice sizes are all set as = 500. The symbols are the numerical data from NEMD simulations. 

IV. CONCLUSIONS 

With this short review, we have briefly introduced the heat conduction of one-dimensional nonlinear lattices and 
the role of renormalized phonons and effective phonon theory. Although there are more than one type of excitation 
modes in the nonlinear lattices, detailed analysis on the temperature dependence of sound velocities reveals that the 
energy carriers responsible for the heat conduction should be the renormalized phonons. 

Within the framework of renormalized phonons, a phenomenological theory, effective phonon theory, has been 
developed to explain the heat transport behaviors for general one-dimensional nonlinear lattices. For lattices without 
on-site potential, the effective phonon theory can quantitatively and consistently predict the scaling exponents of 
temperature-dependent thermal conductivities. In particular, for the FPU-/? lattices, the predictions from effective 
phonon theory are verified by NEMD simulations at both low and high temperature limits. For the special class of Hn 
lattices, the predictions from effective phonon theory are in perfect agreements with the obtained numerical results. 

Finally, we need to point out that our work is closely related to the emerging field of phononics which focuses 
on the manipulation of heat flow and information processing with phonons (161]. For a two-scgmcnt thermal diode, 
the rectification of heat flow with the inversion of a temperature gradient is realized by the temperature-modulated 
overlap of the phonon spectra of two segments. It will be crucial to know the temperature dependence of the phonon 
spectrum in advance. Our analysis of the renormalized phonons due to nonlinearity or temperature enables us to 
provide the key information for the basic design of thermal diodes. On the other hand, our effective phonon theory 
dealing with the analysis of temperature-dependent thermal conductivities can be also very helpful for the design of 
novel phononics devices. The key element of phononics devices is the thermal transistor which relies on an effect of 
negative differential thermal resistance (NDTR). And this NDTR effect is a consequence of the system's very steep 
temperature-dependent thermal conductivity, e.g. n oc T^^ with 77 > 1 at least. Therefore, we hope our work about 
renormalized phonons and effective phonon theory can provide solid theoretical support for the ongoing research filed 
of phononics [la|. 




9 



V. ACKNOWLEDGEMENTS 



Parts of the work were carried out at National Supercomputer Center in Tianjin, and the calculations were performed 
on TianHc-lA. This work has been supported by the startup fund from Tongji University (N. L. and B. L.). 



S. Liu, X. Xu, R. Xie, G. Zhang, and B. Li, Eur. Phys. J. B 85, 337 (2012). 

F. Bonetto, J. L. Lebowitz, and L. Ray-Bellet, in Mathematical Physics 2000, edited by A. Fokas et al. (Imperial College 
Press, London, 2000) , 128 (2000). 

S. Lepri, R. Livi, and A. Politi, Phys. Rep. 377, 1 (2003). 

B. Li, J. Wang, L. Wang, and G. Zhang, CHAOS 15, 015121 (2005). 
J.-S. Wang, J. Wang, and J. T. Lii, Eur. Phys. J. B 62, 381 (2008). 

A. Dhar, Adv. Phys. 57, 457 (2008). 

C. W. Chang, D. Okawa, H. Garcia, A. Majumdar, and A. Zettl, Phys. Rev. Lett. 101, 075903 (2008). 
M. Terrano, M. Peyrard, and G. Casati, Phys. Rev. Lett. 88, 094302 (2002). 

B. Li, L. Wang, and G. Casati, Phys. Rev. Lett. 93, 184301 (2004). 

B. Li, J. Lan, and L. Wang, Phys. Rev. Lett. 95, 104302 (2005). 

C. W. Chang, D. Okawa, A. Majumdar, and A. Zettl, Science 314, 1121 (2006). 

B. Li, L. Wang, and G. Casati, Appl. Phys. Lett. 88, 143501 (2006). 
L. Wang and B. Li, Phys. Rev. Lett. 99, 177208 (2007). 

L. Wang and B. Li, Phys. Rev. Lett. 101, 267203 (2008). 

R. Xie, C. T. Bui, B. Varghese, Q. Zhang, C. H. Sow, B. Li, and J. T. L. Thong, Adv. Funct. Mater. 21, 1602 (2011). 
N. Li, J. Ren, G. Zhang, L. Wang, P. Hanggi, and B. Li, Rev. Mod. Phys. 84, 1045 (2012). 

C. Alabiso, M. Casartelh, and P. Marenzoni, J. Stat. Phys. 79, 451 (1995). 

C. Alabiso and M. Casartelh, J. Phys. A 34, 1223 (2001). 
S. Lepri, Phys. Rev. E 58, 7165 (1998). 

B. Gershgorin, Y. V. Lvov, and D. Cai, Phys. Rev. Lett. 95, 264302 (2005). 
B. Gershgorin, Y. V. Lvov, and D. Cai, Phys. Rev. Lett. 75, 046603 (2007). 
N. Li, P. Tong, and B. Li, EPL 75, 49 (2006). 
N. Li and B. Li, EPL 78, 34001 (2007). 

D. He, S. Buyukdagh, and B. Hu, Phys. Rev. E 78, 061103 (2008). 
J. A. D. Wattis, J. Phys. A 26, 1193 (1993). 

G. Friesecke and J. A. D. Wattis, Commun. Math. Phys. 161, 391 (1994). 
F. Zhang, D. J. Isbister, and D. J. Evans, Phys. Rev. E 61, 3541 (2000). 
F. Zhang, D. J. Isbister, and D. J. Evans, Phys. Rev. E 64, 021102 (2001). 
S. Flach and C. R. WiUis, Phys. Rep. 295, 181 (1998). 

S. Flach and A. Gorbach, Phys. Rep. 467, 1 (2008). 

B. Hu, B. Li, and H. Zhao, Phys. Rev. E 61, 3828 (2000). 

K. Aoki and D. Kusnezov, Phys. Rev. Lett. 86, 4029 (2001). 

H. Zhao, Z. Wen, Y. Zhang, and D. Zheng, Phys. Rev. Lett. 94, 025507 (2005). 
N. Li, B. Li, and S. Flach, Phys. Rev. Lett. 105, 054102 (2010). 

K. Aoki and D. Kusnezov, Phys. Lett. A 265, 250 (2000). 

B. Hu and L. Yang, CHAOS 15, 015119 (2005). 

R. Lefevere and A. Schenkel, J. Stat. Mech.: Theory Exp. , L02001 (2006). 

N. Li and B. Li, Phys. Rev. E 76, 011108 (2007). 

L. Nicolin and D. Segal, Phys. Rev. E 81, 040102(R) (2010). 

S. Flach, M. Ivanchenko, and N. Li, Pramana J. Phys. 77, 1007 (2011). 

N. Li and B. Li, J. Phys. Soc. Jap. 78, 044001 (2009). 

H. Zhao, Phys. Rev. Lett. 96, 140602 (2006). 

F. Piazza and S. Lepri, Phys. Rev. B 79, 094306 (2009). 

C. Giardina, R. Livi, A. Politi, and M. Vassalh, Phys. Rev. Lett. 84, 2144 (2000). 
O. V. Gendelman and A. V. Savin, Phys. Rev. Lett. 84, 2381 (2000). 

Y. Zhong, Y. Zhang, J. Wang, and H. Zhao, Phys. Rev. E 85, 060102(R) (2012). 

G. R. Lee-DadsweU, E. Turner, J. Ettinger, and M. Moy, Phys. Rev. E 82, 061118 (2010). 
B. Hu, B. Li, and H. Zhao, Phys. Rev. E 57, 2992 (1998). 



